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Abstract. We review research on a number of situations where a quantum impurity 
or a physical boundary has an interesting effect on entanglement entropy. Our focus 
is mainly on impurity entanglement as it occurs in one dimensional systems with a 
single impurity or a boundary, in particular quantum spin models, but generalizations 
to higher dimensions arc also reviewed. Recent advances in the understanding of 
^ , impurity entanglement as it occurs in the spin-boson and Kondo impurity models arc 

discussed along with the influence of boundaries. Particular attention is paid to 1 + 1 
dimensional models where analytical results can be obtained for the case of conformally 
invariant boundary conditions and a connection to topological entanglement entropy 
is made. New results for the entanglement in systems with mixed boundary conditions 
5o ■ are presented. Analytical results for the entanglement entropy obtained from Fermi 

liquid theory are also discussed as well as several different recent definitions of the 
\Q • impurity contribution to the entanglement entropy. 
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1. Introduction 

The definition of entanglement entropy is based on dividing space into two regions. 
In many cases the systems under study are homogenuous and this division is purely 
fictitious. However, there has also been considerable activity on studying entangelement 
in inhomogeneous systems, the simplest of which contain physical boundaries or a single 
impurity. There are currently several useful measures of entanglement, here we shall use 
the von Neumann entanglement entropy as defined by dividing a bipartite system in a 
pure state at T = into 2 regions, A and B. From the ground state pure density matrix, 
region B is traced over to define the reduced density matrix pa- In most cases we shall 
take A to include the impurity/boundary. From this the von Neumann entanglement 
entropy [H[2], 



is obtained for a subsystem of size r inside a larger system of size R. There are several 
motivations for this work. 

One motivation is to study models of a qubit interacting with a decohering 
enviroment [3j EJ El E] . Such a system is often represented by a 2- level system, or spin- 
1/2, interacting with an otherwise homogeneous, and often one- dimensional medium 
with gapless excitations. Some versions of this model are equivalent to the Kondo model, 
motivating studies of ground state entanglement of an impurity spin with the conduction 
electrons in Kondo models. This entanglement entropy can be easily expressed exactly in 
terms of the impurity magnetization, which, for many models, has been well-understood 
many years ago [5j [7]. This single site impurity entanglement, which we denote by S; mp , 
is reviewed in section 

Another motivation comes from the thermodynamic impurity entropy, Ing, which 
was calculated by Bethe ansatz [HE] for multi-channel Kondo models and then discussed 
more generally from the viewpoint of Conformal Field Theory (CFT) [10]. General 
quantum impurity models, such as occur in condensed matter physics, were argued to 
renormalize to conformally invariant boundary conditions, and Ing was argued to be a 
universal quantity depending only on the boundary condition. Calabrese and Cardy (C 
& C) [H] argued that, for a CFT defined on the semi-infinite line with a conformally 
invariant boundary condition (CIBC) at the end, the entanglement of a region of length 
r with the rest is given by [12j [131 HH E] : 



lag depends on the CIBC, establishing a surprising connection between thermodynamic 
and entanglement entropy. Here a is a cut-off length scale and Si is a non-universal 
number. Both are independent of the CIBC. In the numerical work one-dimensional 
tight binding models are generally considered in which case a is the lattice constant. 
This raised the possibility that this impurity part of the entanglement entropy, S irap , 
might exhibit universal renormalization group (RG) behavior and this was confirmed 
in studies of spin chains with a boundary magnetic field [T3] and of the Kondo 
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model [HUE]. The connection between the thermodynamic impurity entropy and the 
impurity entanglement entropy is outlined in section [3j 

There is a deep connection between (1+1) dimensional CFT and topological phases 
of gapped (2+1) dimensional systems such as occur in the fractional quantum Hall 
effect. While the entanglement entropy of a region in a gapped 2 dimensional (2D) 
system is expected to grow with the length of its perimeter, it was shown that there 
is an additional universal, length independent " topological entropy" [TTJ [18] , — In T> 
which can be extracted from the corresponding (1+1) dimensional CFT. The connection 
between the boundary entropy of C & C [11] and the topological entropy can be clarified 
by considering a (2+1) dimensional system with boundaries, a "Hall bar", containing 
a point contact. There are gapless degrees of freedom living on the edge of the Hall 
bar described by a CFT. (See Fig. [9j) The point contact may renormalize to a CIBC, 
corresponding to breaking the Hall bar into two pieces and the change in topological 
entropy due to this renormalization can be related to the change in the boundary entropy, 
In g. As we will show, this is defined in terms of the entanglement of a section of the 
Hall bar of length r, containing the point contact with the rest of the (infinite) Hall bar. 
The connection with topological phases is discussed in section 13.31 

An alternative type of impurity-related entanglement entropy has also been studied 
for a ID wire with a point defect, which if relevant, effective breaks the system in two 
at low energies pJH EDI EI] . Rather than studying the entanglement of a finite region 
surrounding the point contact and extracting the In g term in Eq. (II. 2p instead the 
entanglement of one side of the point contact with the other was studied. In cases 
where the point defect is relevant, it was found that this entanglement tends to vanish 
with increasing system size, again verifying that entanglement entropy exhibits RG flow 
behavior. This type of impurity entanglement is discussed in section HI 

A surprising result of numerical studies of entanglement entropy in ID 
antiferromagnets with boundaries was the presence of an alternating term decaying 
away from the boundaries [22l [7]. Although a theory of this is still lacking, it was 
shown numerically to track closely the energy density as a function of distance from the 
boundary and a heuristic understanding was obtained in terms of a local dimerization 
induced by the boundary, related to "resonating valence bonds" . This boundary induced 
alternation in the entanglement entropy is reviewd in section [5j 

Although probably less useful as a model of a decohering environment and not 
related to CFT and universal RG concepts, impurity entanglement entropy has also 
been studied for gapped (1+1) dimensional systems, including dimerized and Haldane 
gap spin chains. This aspect is reviewed in section El 

As a final motivation, we note that recent models for qubit teleportation and 
quantum state transfer using quantum spin chains [23j EU E3 ESI E7J ES EHl [301 [31] 
employ models closely related to the quantum spin models reviewed and in many cases 
rely on properties of the entanglement arising from the impurities reviewed here. 
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The simplest definition of the impurity entanglement is to consider the (single site) 
impurity as sub-system A and the rest of the system as sub-system B. A measure of the 
impurity entanglement is then simply given by the von Neumann entanglement entropy 
of the reduced density matrix for A inside a system of total size R. Since A describes just 
a single site one often refers to this as the single site impurity entanglement. We shall 
denote this quantity by s imp to distinguish it from S imp (r, R) defined in later sections 
by the difference in the uniform part of the von Neumann entanglement entropy for a 
sub-system of extent r with and without the impurity present. Since s imp is concerned 
with a single site such a definition through a subtraction is not possible and the explicit 
r dependence through the size of the sub-system A is absent. 

The single site impurity entanglement, Sj mp , has been studied mainly in 2 different 
settings: The spin-boson model [32} SI El [6] and the closely related Kondo model [7j. 
For the case where the impurity is a s = 1/2 system (qubit) it is easy to see J32J HI [7] 
that: 

* imp 

- 5^(1/2 ±m imp )ln[(l/2± 

± 

"^imp being the magnetization of the impurity in the ground-state. For a system with a 
singlet ground state m imp = and s imp is maximal [3] , the qubit is maximally entangled 
with the rest of the system. For a system with a doublet ground-state (R odd) the 
behavior of m; mp is more interesting and exhibits the usual cross over associated with 
Kondo physics. In this case mj mp was studied, for the usual fermion Kondo model, in 
[33] and [3U [35] for example from which s imp can be derived. 
The spin-boson model is defined by: 

H S b = --^ a * + \°* + H osc + ~cr z 53 \( a i + a J)' ( 2 - 2 ) 

g 

where a x and a z are Pauli matrices and A is the tunneling amplitude between the states 
with a z = ±1. H osc is the Hamiltonian of an infinite number of harmonic oscillators 
with frequencies {uj q }, which couple to the spin degree via {X q }- The heat bath is 
characterized by its spectral function J{uj) = ^^2 q ^q^(^q — oj) = 2nauj, w C w c 
(Ohmic heat bath). Efficient NRG calculations can be performed on this model through 
a mapping [36J to the anisotropic Kondo model. This allowed for rather detailed NRG 
studies [32J of s imp as a function of a. Exploiting known exact results for (cr 2 ) and (a x ) 
in the spin-boson problem it has been shown [1] that S; mp in several limits is a universal 
function of h/T^, where Tk is the Kondo scale. For Tk «C h A it was found [4] that: 

lim s imp (a, A, h) = k 2 (a) ( J In ( A ) (2.3) 
r K «ii«A V n J \± K I 
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with k 2 {a) a known cut-off independent funtion so that Sj mp in this limit is a universal 
function of h/Tx- On the other hand, for h <C Tr <C A it can be shown that: 

lim s imp (a, A,h) = s imp (a, A, 0) - fci(a) f ^- J , (2.4) 

with fci(a) a known cut-off independent function. However, in this case Sj mp (a, A, 0) is 
in general non-universal and only the second term exhibits scaling. The general results 




0.2 A 0.4 0.6 0.8 1 

a 



Figure 1. Si mp (a, A = 0.01w c , h) (E on the figure label) as a function of a for a 
range of h/u> c . Reproduced with permission from [4]. 

are illustrated in Fig. CD where Si mp (a, A = 0.01u; c , h) (E on the figure label) is plotted 
versus a for a range of h/u c . For h = s imp (a) is a monotonically increasing function 
where as for non-zero h a maximum associated with the crossover h ~ Tk occurs. The 
case of a sub-Ohmic heat bath has also been studied [5] . 

The usual Kondo Hamiltonian [371 EH] contains a Heisenberg interaction between 
a s — 1/2 impurity spin, S 1 , and otherwise non-interacting electrons. A simple model 
takes a free electron dispersion relation and a 5 -function Kondo interaction: 

H = J G?r[^{-V 2 /2m)i) + J K S 3 (r)ip j (a/2)tP ■ S}. (2.5) 

(Actually, an ultra-violet cut-off of the 5-function interaction is necessary for the model 
to be completely well-defined.) In the ground-state of the Kondo model the impurity 
spin is screened by the conduction electrons through the formation of a singlet. This 
phenomenon is expected to take place on a length scale: 

tK = v/T K <xe 1 K vJ *\ (2 6) 

where v is the density of states per spin band , Tk is the Kondo scale and v the 
velocity of the fermions. Due to the 5— function form of the interaction Eq. (I2.5P can 
be reduced to a one-dimensional model which can be represented by a lattice model 
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of finite extent R (including the impurity), suitable for numerical studies. As outlined 
above and discussed in more detail in section 13. 2[ we expect the case of R odd to reflect 
Kondo physics and possibly scaling with R/^k- However, it has been shown [7] that Si mp 
exhibits weak scaling violations from the expected R/£k scaling. In particular, m 
was shown to take the form: 

-i 2 



imp 



"limp ~ 



J 



K 



7TV 



\n(R/a) 



(2.7) 



where a is a short distance cut off. Hence weak scaling violations are present in mj mp and 
therefore also in s; mp even though s imp clearly displays the expected crossover related to 
Kondo physics. The presence of such scaling violations is illustrated in Fig. [2j 




Figure 2. (a) Weak scaling violations for the single site impurity entanglement 
entropy Si mp (a) [Eq. (|2.1|) ] and the local impurity magnetization mi mp = (Sf) (b). 
All data are for odd length chains between R = 19 . . . 101. Reprinted from [7j. 



Finally, recent work has focused on the entanglement between the 2 impurity spins 
in the two-impurity Kondo model [HH El HO] as well as the pair-wise entanglement 
between end spins in open ended Heisenberg spin chains [41] . 



3. Impurity entropies from the CFT perspective 

In this section we review the concept of thermodynamic impurity entropy and its 
connection with CFT. We also connect it with topological entropy in 2D topological 
insulators and discuss various applications of these general ideas to specific models: 
spin chains with boundary fields, Kondo model and point contact in a Hall bar. 
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3.1. Thermodynamic Impurity Entropy 

Thermodynamic impurity entropy can be measured experimentally from impurity 
specific heat using the thermodynamic identity: 

C = T§. (3.1) 

The impurity contribution to the specific heat, and hence the entropy, can be measured 
by subtracting off the same quantity for the pure system. Such measurements and 
related theoretical calculations have been performed for many years in Kondo impurity 
systems. 

C (T) for a metal with a dilute random array of magnetic impurities is measured and 
C(T) for the pure system is subtracted off. Dividing by the number of impurities and 
extrapolating to zero impurity density gives the contribution to the specific heat of a 
single impurity, Ci mp (T). From integrating Ci mp (T)/T one can, at least in principle, 
obtain S™, L(T). We consider the infinite volume limit, so that (T) becomes a 
function of T only. We may formally define S^ p (T) by calculations on a ID system of 
length R with a single impurity: 

Sg,(T) = lim [S(T, R) - S Q (T, R)]. (3.2) 



R 



Here S(T,R) is the thermodynamic entropy with the impurity present and Sq(T,R) 
is the same quantity without the impurity. S^ p (T) exhibits interesting T-dependence 
which reflects the RG flow of the Kondo model. It can be calculated with high precision 
from the Bethe ansatz solution of the Kondo model first derived by Andrei [42] and 
Wiegmann [43] . In the limit of a weak bare Kondo coupling (T) w In 2 at T ^> T% 
and S^ p (T) ^ at T < T#. This reflects the fact that the bare coupling of the 
magnetic impurity to the conduction electrons is very weak so that we obtain essentially 
the full entropy of a free spin- 1/2, In 2 at T ^> Tk- However, as the temperature is 
lowered the spin is "screened" i.e. it goes into a singlet state and the impurity entropy 
is accordingly lost. The asymptotic values of S^ p (T) at high and low temperatures are 
characteristic of the RG fixed points of the Kondo Hamiltonian. 
The /c-channel Kondo model 

H = I d 3 r ^[^(-V 2 /2m)^ + J K 5\r)^\a /2)^ • S}. (3.3) 
J i=i 

also has S^ p (T) « In 2 at T ^> Tk- However, it was found, from the Bethe ansatz 
solution, that at T — > it has the limiting value \ng where |10j : 

g = 2cos[7r/(2 + k)) < 2. (3.4) 

Heuristically, g represents a "fractional ground state degeneracy" characterizing the 
non- Fermi liquid ground state of the overscreened Kondo model. 

The interesting behavior of the impurity entropy in the multi-channel Kondo model 
was later shown to be a special case of a general phenomenon in quantum impurity 
models. Many models of this type have low energy descriptions in terms of ID 



Entanglement entropy in quantum impurity systems and systems with boundaries 8 



CFT's. Quite general boundary conditions and boundary interactions are expected 
to renormalize, at low energies, to CIBC's. Cardy showed that generally CIBC's can 
be associated with boundary states, \A>. The conformally invariant partition function 
defined on a strip of length R, at inverse temperature j3, with CIBC's A and B at the 
two ends can be written: 

Z AB = tve- pH ^, (3.5) 

where H AB is the Hamiltonian on an interval of length R with BC's A and B at the two 
ends. Alternatively, we may switch space and imaginary time directions and write: 

Z AB =< A\e~ RH p\B > . (3.6) 

Now the system propagates for time R, under the action of the Hamiltonian defined on a 
periodic interval of length f3 with initial and final states \A > and \B >. The boundary 
states can be expanded in a complete basis of Ishibashi states, associated with each 
conformal tower, a: 

\A >= \a >< aO\A > . (3.7) 

a 

(The Ishibashi states \a > are sums over all descendants with equal weight in left and 
right-moving sectors.) Thus, we may write: 

Z AB = ^2 < A \ a ® >< a ®\ B >< a\e~ RH p\a > 

a 



Z^ 



< A\aO >< aO\B > Xa(e" 47rK/ ^) (3.f 



a 

,th 



where Xa is the character of the a conformal tower. Written in this form it is 
straightforward to extract the impurity entropy. Taking the limit R ^> (3 only the 
highest weight state in the conformal tower of the identity operator contributes, giving: 

Zab -> e nRc/6p < A\00 >< 00\B > . (3.9) 

(Here c is the central charge associated with the bulk CFT.) From this expression we 
obtain the entropy: 

d rm , „ , ttRcT 

Sab — " 

where 



Tb = ^[TlnZ AB ] = ^ + \ng AB (3.10) 



g A B = 9A-gB=<A\00><00\B> . (3.11) 

This consists of the bulk term, independent of the boundary conditions and proportional 
to the system size, R, as well as the boundary term In g which is a sum of contributions 
from each boundary. Using the known boundary state corresponding to the non-Fermi 
liquid ground state of the multi-channel Kondo model we can reobtain the Bethe ansatz 
result for S^ (T = 0) from this general CFT formula. 

In g, the T = thermodymamic entropy, is a universal property of fixed points of 
the boundary RG. It has the interesting property that it always decreases under RG flow 
from an unstable to stable fixed point [TO]. Eq. ( 13. provides an example of this: the 
impurity entropy is In 2 at the unstable fixed point, T 3> T#, and always has a smaller 
value at the stable fixed point which occurs at T = 0. 
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3.2. Boundary term in the entanglement entropy 

We now consider a semi-infinite CFT (r > 0) with CIBC, of type A at r = 0. We 
consider the ground state entanglement entropy for the region, < r' < r, ^(r). We 
might now expect some additional term in S(r), ca, which depends on the CIBC A, but 
not on the length of the region under consideration, r: 

W = Hd +e * + T (312) 

We can argue that ca = hig^, the thermodynamic impurity entropy, by the device of 
considering the entanglement entropy for this system at finite temperature. C&C [TTj 
showed that the generalization of S^(r) to a finite inverse temperature, (3, is given by a 
standard conformal transformation: 



o 



(3 . , f2nr 
sinn 



ca + sx/2. (3.13) 



Tia \ (3 

See also [2]. S^r, /?) is defined by beginning with the Gibbs density matrix for the 
entire system, e~^ H and then again tracing out the region r' > r. Now consider the 
high temperature, long length limit, (3 <C r: 

S ^w + H^) +CA+S i +O{e ' , " ,0y (3 ' 14) 

The first term is the extensive term (proportional to r) in the thermodynamic entropy 
for the region, < r' < r. The reason that we recover the thermodynamic entropy 
when r 3> (3 is because, in this limit, we may regard the region r' > r as an "additional 
reservoir" for the region < r' < r. That is, the thermal density matrix can be defined 
by integrating out degrees of freedom in a thermal reservoir, which is weakly coupled 
to the entire system. On the other hand, the region r' > r is quite strongly coupled to 
the region r' < r. Although this coupling is quite strong, it only occurs at one point, 
r. When r ^> f3, this coupling only weakly perturbs the density matrix for the region 
r' < r. Only low energy states, with energies of order 1/r and a neglegible fraction of 
the higher energy states (those localized near r' = r) are affected by the coupling to the 
region r' > r. The thermal entropy for the system, with the boundary at r = in the 
limit r ^> (3 is: 

SV 1 — > — - — \-\ngA + constant, (3.15) 
Qp 

with corrections that are exponentially small in r/ (3. The only dependence on the CIBC, 
in this limit, is through the constant term, In qa, the impurity entropy. Thus it is natural 
to identify the BC dependent term in the entanglement entropy with the BC dependent 
term in the thermodynamic entropy: 

c A = lng A . (3.16) 

This follows since, in the limit, r > /], we don't expect the coupling to the region 
r' > r to affect the entanglement entropy associated with the boundary r' = 0, ca- 
Note that the entanglement entropy, Eq. (13.141) . contains an additional large term not 
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present in the thermal entropy. We may ascribe this term to a residual effect of the 
strong coupling to the region r' > r on the reduced density matrix. However, this extra 
term does not depend on the CIBC as we would expect in the limit r ^> (3 in which 
the "additional reservoir" is far from the boundary. Now passing to the opposite limit 
(3 — > oo, we obtain the remarkable result that the only term in the (zero temperature) 
entanglement entropy depending on the BC is precisely the impurity entropy, In g A , Eq. 
( II. 2p . We may also consider a finite system, of length R, with ICBC at both ends. In 
the case, where both BC's are the same, A, we expect the generalization of Eq. ( 11.21) 
which follows from a conformal transformation: 



S A (r,R) = ^ln 
o 



2R . /71T\ 
— Sill ( ■ 

7ia 



\RJ 

with the equivalent result for periodic boundary conditions 



+ \ng A + S j-, (3.17) 



S™ c (r,R) = -\n 



R . 


(irr\ 


— sin 




ixa 


K~RJ 



+ Sl . (3.18) 

As usual, a denotes the lattice spacing. A result has not been obtained, as far as we 
know, for the case where the CIBC's are different at the two ends. 

As discussed in the previous sub-section, the thermodynamic impurity entropy, In g, 
is a universal quantity with interesting RG behavior. It is then natural to expect that 
the impurity term in the entanglement entropy will behave similarly. In particular, 
consider beginning with a CFT on the semi-infinite line with CIBC A and then adding 
a small relevant boundary interaction: 

H = H CF t - A0(O) (3.19) 

where <p has an RG scaling dimension y < 1. (y — 1 is the marginal dimension for 
boundary interactions, since the action contains only a time-integral over 0(r, r = 0) and 
no spatial integral.) We expect that, under the renormalization group, the Hamiltonian 
will flow to an infrared stable fixed point, at long length scales, characterized by some 
other CIBC, B. The "g-theorem" [TO, HI] implies that the thermodynamic impurity 
entropy at the stable fixed point obeys gs < gA- Typically the flow between fixed 
points can be controlled by introducing a finite length scale which acts as an infrared 
cut off. This is often done by putting the system in a finite box, of size R. It is 
not obvious what will happen when we have no physical box but introduce a length 
scale I by the definition of the entanglement entropy. Will the impurity part of the 
entanglement entropy exhibit a cross over from g A to gs as we increase /? Will this 
cross over be universal? These questions have been investigated numerically, using the 
DMRG method, in a couple of models. 

The first model considered was the ID transverse field Ising model with a 
longitudinal boundary magnetic field [T5] : 



H = - + (1/2)5?) + M5- (3-20) 

j=0 

The bulk transverse field has been tuned to its critical value, 1/2. A weak longitudinal 
field, hb is applied at the boundary, j = 0, only. This boundary field is relevant, with 
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Figure 3. S(r, R) versus r (R = 500) for the transverse field Ising model calculated 
using DMRG with m = 128. Results are shown for free-free and fixed-fixed boundary 
conditions along with four different values of hj, = 0.01,0.03,0.10,0.30. For the 
numerical DMRG data only every 4th point is shown for clarity. The solid and dashed 
lines represent fits of the DMRG results for free-free and fixed-fixed BC's, respectively, 
to the analytical result Eq. (|3.17j) . For the latter case calculations were performed 
formally with = oo. 



dimension 1/2 and therefore induces a boundary RG flow between the only two boundary 
fixed points in this model, corresponding to free or fixed BC. The values of g are g = 1 
(free) and g = l/y/2 (fixed). DMRG results on chains of length up to 800, keeping 
140 block states, showed quite convincingly that the entanglement entropy crosses over 
from (1/12) ln(r/a) to (1/12) ln(r/a) — (1/2) ln2 as r is increased from small values to 
values larger than a cross-over scale, £, determined by hb as £ oc |/ib| -1 / 2 . In figure [3] 
we illustrate these results by calculations on systems with R = 500. In figure H] we 
also show the entanglement entropy with different BC's at the ends of a finite chain, 
fixed-free and fixed up-fixed down. As far as we know, no analytic formulas have been 
derived for these cases. We note that the expression: 



S(r, R) = ^ 



R . /7rr\ 
— sin — 



In 



/ 7ir\ 
cos — 
\4RJ 



+ constant (3.21^ 



a V R ) 

seems to fit the data quite well in the fixed-free case. 

The second model in which this crossover was studied is the Kondo model [161 U] 
For this model, we may consider the entanglement of a region inside a sphere of size 
r surrounding the impurity, with the rest of space which may either be infinite or 
confined to a larger sphere of size R. The impurity entanglement entropy is defined as 
the increase in entanglement arising when the impurity is added to the system. This 
mimics the definition of the impurity thermodynamic entropy which has been well- 
studied experimentally and theoretically for Kondo systems. Due the (5-function form 




Figure 4. S(r, R) versus r (R = 500) for the transverse field fsing model calculated 
using DMRG with m = 128. Results are shown for free-free and fixed-fixed boundary 
conditions along with free-fixed and fixed (up) - fixed (down) boundary conditions. 
For the numerical DMRG data only every 4th point is shown for clarity. The solid 
and dashed lines represents fits of the DMRG data for free-free and fixed-fixed BC's, 
respectively, to the analytical result Eq. (|3.17p . The calculations with fixed boundaries 
are performed by formally setting hf, = ±oo. 



of the Kondo-interaction, the 3D model is equivalent to a ID model. To see this, we 
expand the electron fields in spherical harmonics. Only the s-wave interacts with the 
impurity spin, the other harmonics being completely free. The entanglement entropy 
may be written as a sum of contributions from each spherical harmonic: 

S = J2 S i,™. (3.22) 

l,m 

Only the s-wave part, So,o is affected by the Kondo interaction so that the total impurity 
entanglement entropy is given by the change in So,o when the impurity spin is added. 
Assuming that the Kondo coupling is weak, as is usually the case in experiments, we may 
integrate out Fourier components of the s-wave electron fields except for a narrow band 
around the Fermi sphere. Linearizing the dispersion relation near the Fermi energy: 

k 2 /2m Rj k F /2m + vp(k — kp), (3.23) 

the model becomes equivalent to a relativistic Dirac fermion defined on the half-line, 
r > interacting with the impurity spin at r = 0: 

H 1D « (iv F /2n) / dr[^ L {d/dr)^ L - ^ R (d/dr)M 
Jo 

+ v F \ K ^ L (0)(a/2)i, L (0) ■ S. (3.24) 
Here Xk oc J and a boundary condition: 

MO) = -MO) (3-25) 
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is imposed on the left and right movers. 

To obtain a model amenable to DMRG studies, we could consider a ID tight-binding 
version of this ID continuum model. However, considerable numerical speed-up can be 
obtained by considering a "spin-only" version of the model. This is based on spin-charge 
separation for ID interacting fermion systems, which follows from bosonization. We find 
that only the spin degrees of freedom of the ID electrons interact with the impurity, 
when it is at the end of the chain. [Things are more complicated when it is not at the 
end. The simplifications at the end arise from the BC of Eq. ( 13.25D .] The spin part 
of the Hamiltonian, the only part involving the Kondo interaction, can be written as a 
perturbed S77(2)i Wess-Zumino-Witten non-linear a-model with Hamiltonian: 

H s = (v F /2ir) [ dr(l/3)[J L ■ J L + J R ■ J R ] + v F X K J L (0) ■ S. (3.26) 
Jo 

Here Jl/r{t) are the spin density operators for left and right movers, with the BC 

Jl(0) = 4(0). (3.27) 

This implies that we may regard Jr{j) as the analytic continuation of Ji,(r) to the 
negative r axis: 

Mr) = J L (-r) (3.28) 
and write the theory in terms of left movers only defined on the interval —R < r < R: 

H s = (v/6n) / drJ L ■ J L + v F X K J L (0) ■ S. (3.29) 

J-R 

Now consider Heisenberg antiferromagnetic S=l/2 chain with one weak link at the 

end: 

R-l 

H = J'^ -S 2 + J2^r- S r+ i. (3.30) 

r=2 

For J' K 1 essentially the same low energy continuum limit field theory, Eq. (I3.26f) is 
obtained except that the Fermi velocity, v F is replaced by the spin- velocity, which we 
call v. This model is considerably more efficient to study with DMRG since there are 
only 2 states per site rather than 4. Actually, a drawback of this model is that there is 
an important marginally irrelevant bulk interaction in the low energy Hamiltonian: 

5H = -(gv/27i)J L ■ J R . (3.31) 

with the positive dimensionless coupling constant g of 0(1). This leads to 
logarithmically varying corrections to all quantities which greatly hinders numerical 
work. To circumvent this problem, it is advantageous to add a second neighbor 
interaction, considering instead the Hamiltonian: 

R-l R-2 

H = J' K (S\ ■ §2 + J2S1 ■ Ss^j + S r ■ S r+ i + J2 S r ■ S r+ 2- (3.32) 

r=2 r=2 

For J 2 > J2 ~ 0.2412 the model goes into a gapped dimerized phase, of which the 
exactly solvable Majumdar-Ghosh model with J2 = 1/2 is a special simple case. The 
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gap is driven by the marginal coupling constant g which changes sign at J°, becoming 
marginally relevant. For J 2 < J 2 t ne model remains gapless with the marginal coupling 
constant g and the spin velocity v varying smoothly. Right at the critical point, J 2 = Jf , 
g = 0. At this point all logarithmic corrections vanish and it becomes possible to 
extract meaningful results from numerical studies of relatively short chains. Therefore, 
we largely focussed on this J 2 = J| model. As discussed in more detail in Ref. [15] and 
references therein, we then see that the low energy effective field theory description of 
the spin only model, Eq. (13.321) . with J 2 = J| (g = 0) is the same as that of the usual 
electronic version of the Kondo model. 

As summarized in sub-section 13.11 . the thermodynamic impurity entropy, S^ p (T) 
decreases monotonically from m In 2 at T ^> T# to zero at T Tk- We might expect 
that the impurity entanglement entropy would behave the same same with the energy 
scale, T replaced by v/r where r is the size of the region A. While we ultimately 
confirmed this result two interesting subtleties were encountered en route. 

First of all, even for an infinite system size, R, we found that the entanglement 
entropy has an alternating term, Sa, which decays only slowly with r |22j : 

S(r, R) = S V (J' K , r, R) + (-l) r S A (J' K , r, R), (3.33) 

We discuss this in section El Although we expect Sy to be essentially the same in the 
spin chain Kondo model as in the fermion Kondo model, the same is not true of S A - 
(More generally, the fermion model is expected to have a term in the entanglement 
entropy oscillating at wave-vector 2k F .) Henceforth, in this section, we focus on the 
uniform part, Su only. 
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Figure 5. DMRG results using spin-inversion with m = 256for 5 l i mp (0, r, R) for 
TZ even (R = 100,400). Results are shown for the J\ — J 2 spin chain model at J|. 
Inset: DMRG results for 5i mp (l, r, R) for R odd and even. Data for R = 100, 101 and 
R = 400, 401 are indistinguishable. Si mp (l, r, R) appears to vanish in the scaling limit 
r — > 00 with r/R held fixed. 
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We extracted an impurity part from Su by: 

S imp (J' K ,r,R) = Su(J' K ,r,R) - Su(l,r - 1,R- 1), r > 1. (3.34) 

Note that we are subtracting the entanglement entropy of a chain where all couplings 
have unit strength and 1 site is removed. This corresponds to subtracting the 
entanglement entropy of the system without the impurity. It is important here that 
we do not subtract the entropy for the same values of r and R but with J' K = since, 
as we discuss below, entanglement with the impurity can survive, even in this limit, 
when R is even, in the spin-singlet ground state. (See Fig. [5]). This is related to the 
second subtlety that we encountered: a very strong dependence of S imp , as defined by 
Eq. (13.341) . on the parity of R, even after subtracting off the part alternating in r. This 
is illustrated by some of our DMRG data shown in Fig. El 




Figure 6. Universal scaling plot of S- lmp for fixed r/R, (a) for R, < 102 even, (b) for 
R < 101 odd. DMRG results with m = 256 for the J\ — J2 chain at Jf for various 
couplings J' K . The lines marked 7r^/(12r) are the FLT prediction: Eq. (|3.39p . (c): the 
location of the maximum, (r/^) max , of S{ mp for odd R, plotted versus r/R. Reprinted 
from [7j. 

This figure tests the conjecture that the impurity entanglement entropy shows 
universal RG scaling behavior. We find that the data for the impurity entanglement 
entropy for various Kondo couplings can be collapsed onto scaling curves which depend 
only on the dimensionless ratios r/^K, where is the Kondo length scale, and r/R. 
However, there are two different sets of curves depending on the parity of R. Note that 
the curves differ markedly for but elsewhere look similar. Indeed it looks likely, 

and presumably must be the case, that as r/R — » 0, the curves become identical for 
even and odd R. Focusing on the even R curves, which seem close to the r/R — > limit, 
the scaling curves seem to be approach a monotone decreasing function with the value 
In 2 at r <C and zero for r 3> This is exactly what we would expect from the 
RG theory of the Kondo model and mirrors the T-dependence of the thermodynamic 
impurity entropy, reviewed in subsection (3.1). In particular, \ng = In 2 at the short 
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distance weak coupling fixed point, corresponding to a paramagnetic spin-1/2 impurity 
but In g = at the long distance strong coupling fixed point, corresponding to the spin 
being screened. 

A qualitative understanding of the surprising difference between even and odd R 
can be obtained using a "resonating valence bond" picture of the ground state of the 
s — 1/2 Heisenberg antiferromagnetic chain [7|. Consider first even R. Any singlet 
state can be written as a linear superposition of products of singlet states formed by 
pairs of spins. Conventionally one draws a line or "valence bond" between pairs of spins 
contracted to a singlet, (| fl> — \ IT>)/v / 2- It can easily be shown that by restricting 
to terms in which none of the lines cross we get a complete linearly independent set of 
singlet states for a S=l/2 chain. Furthermore, by adopting a convenient sign convention 
for valence bond states, it can be proven that all terms in the sum have non-negative 
coefficients. We may heuristically associate the impurity entanglement entropy with the 
valence bond originating from the impurity spin at site 1. If this spin forms a singlet 
with a spin at a site inside region A (at a site with index < r) then we consider this 
not to contribute to the impurity entanglement entropy, Si mp of region A. On the other 
hand, if site 1 is paired with a site outside region A (with index > r) then we consider it 
to contribute to S^mp. The ground state is a sum over many valence bond configurations 



Figure 7. Typical impurity valence bond configurations for R even (a) and R odd 
(b). Note the unpaired spin on the 7 th site in (b). Reprinted from [7J. 

so we could imagine relating Si mp to the probability of this "impurity valence bond" 
(IVB) extending out of region A [161 E]- (See Fig. 0) Similar ideas have been explored 
in [35] . As J' K gets smaller, the IVB gets longer ranged. This follows because, in general, 
valence bonds tend to be nearest neighbor, or at least quite short range in order to take 
advantage of the nearest neighbor antiferromagnetic interactions. However, as J' K gets 
weaker there is less and less energetic advantage in a short impurity valence bond. In 
the extreme case J' K = the impurity valence bond extends with significant probability 
over the entire chain. We expect that the typical length of this impurity valence bond 
should be the Kondo screening cloud size. This is a characteristic length scale, of 
order v/Tjc, at which the effective Kondo coupling becomes of order one. In a metal, 
one heuristically imagines an electron in a quasi-bound state forming a singlet with 
the impurity spin with being the extent of the bound state wave-function. In the 
spin chain realization of the Kondo model this quasi bound state corresponds to the 



Entanglement entropy in quantum impurity systems and systems with boundaries 17 

impurity valence bond. Thus the fact that S imp starts to decrease when r exceeds 
is very natural from this Kondo screening cloud viewpoint. Unfortunately, it seems 
extremely difficult to make this more than a heuristic argument. A problem is that 
the valence bond basis, while complete and linearly independent, is not an orthonormal 
basis. Although one could define a probability distribution for the length of the IVB 
it is hard to relate this to physical quantities such as the correlation length or the 
entanglement entropy. However, we can formally write: 



with p the probability that the IVB connects the impurity spin to a spin in region A. 
This follows, since if the IVB does not cross the boundary between regions A and B 
its contribution to S; mp is obviously zero, on the other hand, if the IVB connects to a 
spin in region B it will contribute a factor of ln(2). It is possible to be a little bit more 
quantitative using the recently introduced valence bond entropy [171 HHj . If one simply 
focuses on the IVB connecting the spin impurity with the rest of the system, we expect 
the probability V that the IVB has a length r to decay like ^k/^ 2 in the regime r ^> 
thus giving 



Such a 1/r 2 behavior is expected for a pure system [19] . 

Now consider the case of odd R. The ground state now has a total spin of 1/2. 
We may again represent it by valence bonds but there are (R — l)/2 of them with 
one unpaired spin. This unpaired spin may or may not be the impurity spin. As J' K 
gets weaker it becomes more and more likely that the impurity spin is unpaired. Clearly 
when J' K — the other R—l spins form a singlet leaving the impurity in a paramagnetic 
state. In this limit, for odd R, S imp is precisely zero, due to the way we have defined 
it. Numerically, we find that the maximum in Si mp , for odd R, occurs at a value of 
J' K such that £k oc R. This is due to trade-off between two competing effects. As we 
decrease £k from large values we increase the probability of having an IVB (i.e. of the 
impurity not being the unpaired spin). This tends to increase S imp . However, once £x 
becomes less than R, the shortening of the IVB with decreasing becomes important 
and decreases Simp. It is now clear that the limits oo and R — > oo don't commute. 

Taking R — > oo for any fixed r and £x gives the same S imp (r/£x) as occurs for even 
R: a monotone decreasing function. However, holding R fixed and varying £k gives a 
maximum of Si mp in the vicinity of £x ~ R. Again, this is largely a heuristic picture. 

While we so far only have heuristic descriptions of Si mp when is of order R and/or 
r, CFT methods can be used to calculate an analytic expression for Si mp in the opposite 
limit £k <C r (for any ratio oir/R) [T6l[7j. (Note that in this limit the dependence on the 
parity of R disappears.) This calculation is based on doing perturbation theory for Si mp 
in the leading irrelevant operator at the strong coupling, infrared fixed point. In fact 
such perturbation theory is very powerful and very well-known for the Kondo effect, 
going under the name of Nozieres local Fermi liquid theory (FLT). It has been used 



-IVB 
imp 



(1-p) In 2 



(3.35) 




(3.36) 
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long ago to calculate the leading dependence at low temperature of thermodynamic and 
transport quantities. This approach can be based on the continuum limit field theory 
of Eq. ( 13.261) . The infrared stable, strong coupling fixed point Hamiltonian does not 
contain the impurity spin since it is screened and breaking this singlet costs an energy of 
0(Tk)- The low energy Hamiltonian at energy scales -C Tk only contains the continuum 
WZW fields. In this simple, spin-only model, the only effect of the Kondo interaction, 
once it has renormalized to strong coupling, is to switch the finite size spectrum of 
the WZW model between the s = and s — 1/2 conformal tower, corresponding to 
removing one site. To perturb around this fixed point we must identify the leading 
irrelevant operator which we expect to appear as a boundary interaction at r = only. 
[It must be understood that this low energy theory is only valid at length scales large 
compared to We may think of the boundary interaction as being smeared over a 
distance of 0(£k) but this is effective the same as being right at the boundary, r = 
in the effective Hamiltonian.] The leading irrelevant boundary operator, which must 
have SU{2) symmetry, is simply J\{r = 0). [Recall that the BC Jl(0) = Jr(0) means 
there is only one boundary current operator.] Very fortunately, this leading irrelevant 
interaction is proportional to the bulk energy density, as we see from Eq. (13. 29ft . Defining 
this energy density: 

H s , L (r) = (v/Qtt)J l -J l (3.37) 

the leading irrelevant interaction at the infrared stable fixed point is conventionally 
written: 

H int = -«k)H s , l (0). (3.38) 

This can be taken as a precise definition of the crossover length scale The fact 
that H int can be written in terms of 7i. St L is very convenient for calculating the leading 
perturbation to the entanglement entropy because we may simply take over the earlier 
results of Calabrese and Cardy. These imply the leading correction: 

S imp ^7rW(12r), (3.39) 

representing first order perturbation theory in valid when r 3> £k- This result is 
obtained for infinite R but we may obtain the finite R result by a standard conformal 
transformation: 



imp ~ 12R 



r Tcr 
l + vr(l--)cot(- 



(3.40) 



We find that this formula fits our numerical data, for both even and odd R, extremely 
well as illustrated in Fig. [HI Note that there are essentially no free parameters in the 
fit except for which can be determined independently as a function of J' K and is 
expected to behave exponentially at weak coupling: 

^ocexp[1.38/J^]. (3.41) 

(The parameter, 1.38 is not a fitting parameter but is rather determined from a careful 
mapping of the weak coupling of the end spin to the Kondo coupling in the usual 
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400 



Figure 8. DMRG results with m = 256 for the 3\ — Ji chain. Si mp (J' K = l,r,R) for 
R = 400,401 with J 2 = J| compared to the FLT prediction, Eq. (|5T30|l . 



fermion Kondo model.) This good fit of the CFT predictions to DMRG results on the 
spin chain version of the Kondo model is rather striking confirmation of the universality 
of entanglement entropy. 

3.3. Topological entanglement entropy 

We now consider a class of gapped 2D systems known as topological insulators. These 
occur most famously as models for the fractional quantum Hall effect although other 
experimentally relevant possibilities have been conjectured. The gapped excitations 
of these systems can exhibit non-abelian statistics and are currently of great interest 
as possible topological quantum computers. It is rather difficult to demonstrate 
numerically that a given microscopic model has a topological ground state of a given 
type. Recently, it was observed that this information can be extracted from the 
entanglement entropy. If we consider the entanglement entropy of a finite region inside 
a 2D insulator of infinite extent then we expect "area law" behavior: 

S = aR' (3.42) 

where R' is the length of the perimeter of the region and a is a non-universal constant. 
However, it was recently shown [T71 [18] that there is a sub-leading universal topological 
term in the entanglement entropy which is independent of the length and shape of the 
perimeter: 

S = aR'- In V. (3.43) 

To actually extract this term, it was proposed to divide the infinite 2D space into 3 
or 4 imaginary finite regions and 1 infinite one and to calculate a sum and difference 
of various entanglement entropies so that the term oc R cancels. The connection of 
this term with the impurity entropy, In g was suggested in [IT]. (See also [IB]). It was 
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further elucidated in using the connection of a 2D topological insulator with a ID 
edge model. Consider for example a "quantum Hall bar", a macroscopic sample of a 
topological insulator with edges. It is known that the electric current responsible for 
the Hall conductivity, o~ xy , flows around the edges only, in a clean sample, in a direction 
determined by the magnetic field direction, which is perpendicular to the plane. The low 
energy excitations on the edge correspond to a chiral CFT, meaning that the excitations 
are moving in only one direction. By considering a long thin hall bar, such as in Fig. 
([9]) one can formally group together the right moving excitations on the upper edge and 
the left moving excitations on the lower edge to obtain a parity symmetric CFT. (As 
in all CFT's the left and right moving sectors are decoupled.) This is a particularly 
convenient formulation if a constriction is created in the Hall bar, corresponding to a 
narrowing of the bar (usually imposed with gate voltages) at one point, as shown in 
Fig. (EJ). The constriction leads to "back-scattering" i.e. reflection of right-movers 
approaching the constriction on the upper edge into left-movers leaving it on the lower 
edge. Integrating out the gapped bulk modes, this can be described by a purely ID 
parity symmetric field theory: a CFT with a local back-scattering interaction. In some 
situations, for example a simple CFT corresponding to spinless fermions with repulsive 
interactions, this back-scattering is a relevant interaction and can block all transport 
between left and right sides of the ID system at low energies and long lengths. (This can 
be understood in terms of the boundary RG discussed in subsection 2.1). To make this 
connection more explicit let's consider the case where the constriction has zero width 
in the x-direction and is parity invariant so that it acts only on the parity even channel 
of excitations. It is then possible to make a "folding transformation" mapping the ID 
infinite system to a ID semi-infinite system, r > with the scatterer at r = 0. The 
incoming parity even excitations are mapped to left movers and the outgoing parity 
even excitations to right movers. The back-scattering interaction becomes a boundary 
sine-Gordon interaction for the boson field corresponding to the parity even excitations. 
When the back-scattering is relevant it leads to an RG flow from Neumann to Dirichlet 
BC's on the parity even boson. This is associated with a non-zero change in the impurity 
thermodynamic entropy, In g which depends on the strength of the bulk excitations. In 
the original 2D model, we may think of the relevant constriction as effectively breaking 
the Hall bar into 2 pieces, each with its own chiral edge modes, as shown in Fig. ([9]). 

Again we may relate the impurity thermodynamic entropy to the impurity 
entanglement entropy. To do this, in the original 2D model, it is convenient to 
consider the entanglement of some region, A, containing the constriction and extending 
a distance r on either side of it, with the rest of the infinite Hall bar. The corresponding 
entanglement entropy, in the ID folded formulation can be decomposed into a sum of 
contributions from even and odd parity modes. Only the even parity part is affected 
by the constriction. The corresponding entanglement entropy for the even parity edge 
modes is: 




(3.44) 
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Figure 9. A long thin Hall bar with a chiral edge current and a constriction. In 
sub-section 3.3 we let R — ► oo and consider the entanglement entropy for the region of 
width 2r with the rest of the Hall bar. In subsection 3.4 we take R finite and consider 
the entanglement of the left side of the entire system with the right side. 

the standard result for a finite region in an infinite CFT with c = 1 (eg. with periodic 
BC's). Note that we have so far only considered the entanglement between region A 
and the rest due to the gapless edge modes. We expect additional contributions from 
the bulk modes, 2aw, where w is the width of the Hall bar and hence 2w is the length 
of the boundary of region A. Note that the extra term (1/3) lnr does not occur in the 
formulation of [T71 [18] where the region A does not include any physical boundaries and 
hence does not include any edge states. 

Now consider the change in entanglement entropy for region A with the rest, as 
we increase the length, 2r, of region A. We then expect an RG flow from the clean 
Hall bar at small r to a Hall bar broken in two at large r, with the crossover length 
scale determined by the strength of the back-scattering. From the viewpoint of the low 
energy ID model, the change results from the RG flow of the BC and is given by the 
change in \ng: 

Sir - Suv = hi^iR - ln^uy. (3.45) 

[Note, in particular, that the (1/3) lnr term is the same in either perfectly transmitting 
or perfectly reflecting limits. In the first case it is the bulk ID entanglement entropy for 
a region of length 2r in an infinite system. In the second it is twice the entanglement 
entropy for a region of length r in a semi-infinite system.] On the other hand, from 
the viewpoint of topological entropy, the constriction has broken the Hall bar into two 
pieces, without effecting the length of the perimeter, separating region A from the rest, 
2w. This is expected to double the topological entropy: 

S m -S vy = -\nV, (3.46) 

since, in the infrared limit, we simply get twice the entanglement entropy of the Hall 
bar in the region < r' < r with the region r' > r, 2(aw — \nV) with zero entanglement 
between left and right sides, through the constriction. Thus we conclude that: 

In V = In guv - In g m . (3.47) 

Thus the topological entropy of a gapped 2D system is equal to the change in 
impurity entropy of the corresponding ID gapless edge theory under a change in CIBC 
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corresponding to breaking the system into two pieces. We note in passing that more 
generally, the local interactions associated with the constriction could lead to non-trivial 
boundary conditions not simply corresponding to perfect transmission or reflection. This 
would correspond to a different value of gi^ and probe other features of the topological 
insulator. 

4. Bulk impurity effects 

Another novel type of impurity entanglement entropy was studied in [19] and [20, 12T] . 
We consider the same type of model as discussed in sub-section 13.31 a Luttinger liquid 
with a back-scatterer, considering the ID formulation of the model. We now let the total 
system size 2R be finite with the back-scatterer in the centre. Rather than considering 
the entanglement of a central region containing the origin with the rest, these authors 
considered instead the entanglement entropy of the region to the left of the constriction 
with the region to the right. With no back-scattering this is: 



With relevant backscattering, the ID wire is effectively broken into two disconnected 
parts at long length scales so we might expect the entanglement entropy to vanish 
asymptotically. This was studied numerically using DMRG in a critical Heisenberg 
XXZ spin chain with Hamiltonian: 



For uniform couplings, this model is in a gapless Luttinger liquid phase for — 1 < A < 1. 
Here Jj = 1 for all links except for one in the middle of the chain where it has the value 
<^imp < with < J imp < 1. This weakened link is a relevant perturbation for A > but 
irrelevant for A < 0. The DMRG results were consistent with S/\nR going to zero at 
large R for < A < 1 but going to (1/6), the expected c = 1 value, for —1 < A < 0. A 
quantitative theory of the independence has not yet been developed, as far as we know. 

Properties of the impurity entanglement in the s = 1/2 XX chain in a transverse 
field where the field strength is perturbed at a single central site have also been 
investigated [5U E2] . 

5. The alternating part of the entanglement entropy 

5.1. Open boundary induced alternation 

As discussed in Ref. [531 122] . in the case of AF spin chains with OBC, the von Neumann 
entropy can be written as a sum of two contributions: 




(4.1) 




(4.2) 



j 



S(r,R) = S v (r,R) + (-l) r S A (r,R). 



(5.1) 
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Figure 10. Schematic picture of a chain of total length R with OBC. The system has 
been cut in two parts A and B. Sub-system A. of length r, includes one open end. 

The uniform part Sy(r, R), in the case where the block of size r contains one open end 
(depicted in Fig. [TO]), is given by the CFT result [UJ, Eq. fETTTj) : 
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Figure 11. Results from ED of XX chains with R = 2000 spins 1/2. The von 
Neumann entropy S(r, R = 2000) is plotted vs the subsystem size r for PBC (upper 
symbols) and OBC (lower symbols); the insets shows zooms around the chain center. 
The full lines are fits: Eq. jOI with s x = 0.726 for PBC. In the OBC case, uniform 
and staggered terms °'° 5 r~ ( '~/ 1 } t .\' 25 have been added to Eq. ()5.2j) . Reprinted from [22] . 

where In g is the boundary entropy introduced in Ref. [10] and si is a non-universal 
constant. S(r, R) can be exactly computed numerically using either DMRG or exact 
diagonalization (ED). For the S = 1/2 XXZ chain, 

R-l 

^xxz = J ^ ^ (S r S r+ i + S^S^i + A.S r S r+ ij , (5-3) 

which is critical for |A| < 1, there is a special point at A = (XX point) where the 
spin chain is equivalent to free-fermions. There, one can numerically compute the von 
Neumann entropy [5l] over for very large systems using ED, as shown in Fig. [TT1 both 
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for PBC and OBC with R = 2000 sites. For periodic chains, the numerical results for 
the entropy are very well described by the expression [11], Eq. (13.181) : 



S(r, R) 



In 
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CKT\ 


— sin 




an 


<~r) 



+ Si, 



(5.4) 



with c = 1 and si ~ 0.726 as predicted in Ref [55] for free-fermions. On the 
other hand for OBC, besides the uniform logarithmic increase Eq. (15.21) . there are 
additional uniform and staggered terms which decay away from the boundary (0.055 — 
(-l) r 0.25)/ [f sin (f )]. Not predicted by CFT, the origin of the alternating term has 
been carefully investigated using ED and DMRG for critical XXZ chains in [22]. Such 
a phenomenon has also been observed in several other cases where DMRG were applied 
for open systems. As discussed above, in Ref. [20l |21~] Peschel and co-workers, studying 
the effect of interface defects in critical spin chains, reported the observation of such 
oscillations. Fermions or bosons confined in ID geometries with OBC are also affected by 
such a modulation, as reported for fermionic [2E E3, EB] and bosonic [55] Hubbard-like 
models. 




Figure 12. von Neumann entropy computed by DMRG for the 5* = 1 bilinear- 
biquadratic model Eq. (|5-5[) on open chains for two critical points: (a) 9 = — 7r/4 and 
(b) 9 = 7r/4. Black lines are fits to the expression Eq. (|5.2j) with (a) c = 1.5 and (b) 
c = 2. Reprinted from |56j . 



An interesting example of critical spin chain is the S = 1 bilinear-biquadratic model 



H 



R-l 

E 

r=l 



COS 9 ( S r ■ S, 



r+1 



sin 6 (S r ■ SV+i 



(5.5) 



which displays conformally invariant critical points at 9 = it/ '4 with c = 2 and at 
9 = — 7r/4 with c = 3/2. As studied by Legeza and co-workers [56], the dominant 
correlations at k = 7r for 9 = tt/4 and k = 27r/3 for 9 = —n/4 show up in the 
von Neumann entropy as displayed in Fig. [121 However, a quantitative study of this 
boundary-induced term was not achieved in Ref [56] where the authors simply performed 
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a fit to the expression Eq. f!5.2p restricting to the lower points: r = 2p (0 = —tt/4:) and 
r = 3p (6 = 7r/4), with p integer. Another interesting case has been studied with 
DMRG by Roux and collaborators [HE] in the context of spin 3/2 fermionic cold atom 
with attractive interactions confined in a ID optical lattice. They also found OBC- 
induced non-uniform features in the von Neumann entropy S(r, R) with 2kp oscillations 
that also appear in the local density and kinetic energy t(r, R). Oscillations of S(r, R) 
and t(r, R) appear to be directly related as displayed in Fig. [131 Varying the parameters 
of the fermionic Hubbard model studied in this context [58], there is a critical point with 
c = 3/2 which separates two phases with c = 1. Imposing the non-uniform part of the 
von Neumann entropy to be directly proportional to the non-uniform part of the kinetic 
energy a good fit was obtained by Roux and collaborators with a clear jump in the 
central charge c = 1 — > 3/2 at the critical point [58] . 



2 - 
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32 64 96 128 



Figure 13. Von Neumann block entropy S(x) for a block of size x (a) and local kinetic 
energy t(x) (b) around the critical point of a spin 3/2 fermionic Hubbard model with 
attractive interactions U/t = —1.2 for fixed V/t = —2 and n = 0.75 with L = 128. Fits 
using the oscillations of the local kinetic energy term (b) are quite accurate, allowing 
for the determination of the central charge c. Reprinted from [58] . 

5.2. Valence bond physics 

Open boundary induced oscillations in the entanglement entropy appear to be a quite 
general phenomenon, as also observed in a valence bond physics framework [4TI l4"8l [60| 
1ST] . For SU(2) invariant spin systems, the Stot = sector can be studied by quantum 
Monte Carlo simulations in the valence bond basis [62]. In such a framework, one 
can define a Valence Bond Entropy which displays surprising similarities with the von 
Neumann entropy [4"?1 l4"8 l |60 | I6T]. For random bond spin chains, Refael and Moore [16] 
achieved a very nice calculation of the von Neumann entropy, simply observing that for a 
single valence bond configuration (as depicted for instance in Fig. 03] (a)), the entropy is 
just given by N(r) x In 2, where N(r) is the number of singlet bonds crossing the interface 
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Figure 14. (a) Schematic picture for a particular valence bond state favored by OBC. 
The number of crossing bonds N(r) is given, (b) Corresponding von Neuman entropy. 

between the two sub-systems. Such an idea was successfully checked numerically in the 
random bond case [Ml EH HZj- In the disorder free case, while the ground-state is 
a highly non-trivial superposition of a huge number of valence bond configurations, 
such a phenomenological approach appeared to be extremely useful to understand the 
oscillating features [22]. Indeed, since the translational invariance is explicitly broken 
by the open ends, there is tendency towards dimerization near the open edges. Such 
an effect can be computed very precisely (see below) but already at a qualitative level, 
this short-range singlet formation in the vicinity of open boundaries can be interpreted 
as an alternation of strong and weak bonds along the chain, thus leading to a similar 
alternation of S(r, R). Indeed, the boundary spin at r = 1 will have a strong tendency to 
form a singlet pair with its only partner on the right hand side. On the other hand the 
spin located at r = 2 will be consequently less entangled with its right partner at r = 3 
since it already shares a strong entanglement with its left partner. A typical valence 
bon configuration favored by OBC is depicted in Fig. dH (a) with the corresponding 
entropy (b). Such a qualitative interpretation in term of "weak-strong" modulation is 
directly related to OBC induced Friedel-like oscillations that one can investigate in a 
more quantitative way, as we do now. 

5. 3. Entropy oscillation and dimerization for critical XXZ chains 

The alternating part Sa{t, R) has been studied in Ref. [22] all along the critical regime of 
the XXZ chain |A| < 1. DMRG results for S(r, R) = S v (r, R) + (-l) r S A (r, R) are shown 
in Fig. [15] for R = 100 and various A. One sees immediately that the oscillating part 
varies with A and decays faster in the ferromagnetic regime. More quantitatively, it was 
shown [22] that the alternating part Sa(t, R) is directly proportional to the alternating 
term in the energy density. The energy density for XXZ spin chains: 

(hr) = (S*S* +1 + S y r S y r+1 + AS*S* +1 ) (5.6) 

is uniform in periodic chains. On the other hand, an open end breaks translational 
invariance and there will be a slowly decaying alternating term or " dimerization" in the 
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Figure 15. Von Neuman entropy of XXZ chains of lengths R = 100 with OBC 
computed with DMRG for various anisotropics A. 



energy density 

(h r ) = E v (r) + (-l) r E A (r), (5.7) 

where Ea(t) becomes nonzero near the boundary and decays slowly away from it. E A {r) 
is obtained by abelian bosonization modified by OBC [65J. In the critical region |A| < 1, 
one gets [221 U\ 

E A (r,R)cx 1 (5.8) 

[f sm (l)] 

where K is the Luttinger liquid parameter defined as K = 7r/(2(7r — cos" 1 A)) so that 
K — 1 for an XX spin chain, K = 1/2 for the AF Heisenberg model, and K — ► oo 
for the ferromagnetic Heisenberg case. Based on DMRG data obtained [22] on critical 
open chains of sizes 200 < R < 1000, we find a proportionality between S A and Ea- 
More precisely, plotting Sa as a function of — Ea for various values of the anisotropy 
A in Fig. [161 we find a linear relation Sa = —oEa with a prefactor perfectly described 
by a — \ij sin fi, as shown in the inset of Fig. [161 with \i = cos -1 A. We note that the 
velocity of excitations for the XXZ model is given by v = 7r(sin/x)/(2//) so that we may 
write this relation as 

S A = -(na 2 /2v)E A , (5.9) 

where we have introduced the lattice spacing, a, to make the entanglement entropy a 
dimensionless quantity (Ea has dimensions of energy per unit length.). We emphasize 
that Eq. (15.91) even holds for the Heisenberg model (with a proportionality coefficient 
a = 1) where both E A and Sa display the same logarithmic corrections. Indeed, at the 
Heisenberg point, A = 1, Eq. (15. 8p will have some logarithmic corrections due to the 
presence of a marginally irrelevant coupling constant in the low energy Hamiltonian, 
leading to [22l[7j Ea(t) oc l/[A/r(lnr) 3 / 4 ] in the limit R — > oo. It is highly non-trivial to 
include both the log corrections and the finite size effects in Ea(t, R). However, there is 
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Figure 16. Linear behavior of the alternating part of the entanglement entropy, Sa 
as a function of the alternating energy density. —Ea, both computed using DMRG on 
critical open XXZ chains of size 200 < R < 1000 for a few values of the anisotropy 
A. Data from ED at A = are also shown for R = 2000. Dashed lines arc linear fits 
of the form Sa = —&Ea (see text). The inset (a) is a zoom close to 0, showing data 
for A = —3/4 and —1/2. The inset (b) shows the prefactor a vs A extracted from 
the numerical data (circles), for a larger set of values of A, which is compared with 
ir/2v — ^/sin/i, /i = cos -1 A. Reprinted from [22] . 



a simple result at r = R/2. Including the cubic term in the /3-function for the marginal 
coupling constant, and other higher order corrections [66], this becomes: 




where ao, «i and a2 are constants. If one allows a frustrating second neighbor coupling 
J2 in the chain, at J% = 0.241167 J, this model is at the critical point between gapless 
and gapped spontaneously dimerized phase and the marginal coupling constant, and 
hence the log corrections are expected to vanish here. In both cases (J2 = and 
J 2 = Jg), we found proportionality between Sa and Ea, as shown in Fig. [T71 The sum 
Ea(R/2, R) + aSA(R/2, R) is found to rapidly decay as a power-law, with a power ~ 2.5 
(see Fig. [T71) . Interestingly, for J2 = 0.241167, again linearity is observed, but with a 
prefactor a ~ 1.001689 not related to the spin velocity, v, which we have determined to 
be v ~ 1.174(1) [45]. Hence, Eq. (15. 9p does not hold in this case. 

We emphasize that this alternating term in S(r, R) is universal and should not be 
regarded as a correction due to irrelevant operators. First of all, it is not a "correction" , 
since it is alternating. Secondly, it decays with the same power law as EA{r, R) which is 
seen to be a property of the fixed point, not the irrelevant operators. (However, for the 
Heisenberg model, A = 1, the log factor in EA{r, R) is due to the marginally irrelevant 




0.08 





(5.10) 




Figure 17. (a) Entanglement entropy S(r) (black squares) and alternating part of 
the energy density Ea{t) (open circles) computed with DMRG at Jf for R = 200 
sites. The dashed line is Eq. ()5.8|) . The uniform part of S(r), obtained by taking 
S(r) + (— l) r aEU(?~) with a = 1.001699, is represented by open squares. The best fit, 
shown by a red curve, is indicated on the plot, (b) Comparison between the alternating 
part Sa and Ea [Eq. (|5.9|) ] from DMRG with m = 512 states, for Heisenberg models. 
Power-law decay of SU(i?/2, R) + aEA(R/2, R) drawn in a log-log plot, with a = 1 for 
the nearest neighbor chain (J2 = 0) and a = 1.00169 at the critical second neighbor 
coupling J 2 = 0.241167. Lines are power-law fits: ~ i?" 256 for J 2 = 0.241167 and 
R~ 2 - 59 for J2 = 0. (a) and (b) are respectively reprinted from [7] and [2"2~| . 



operator.) The presence of a universal alternating term in S(r,R) is connected with 
the antiferromagnetic nature of the Hamiltonian (not appearing, for example, in the 
quantum Ising chain [15]) and does not seem to follow from the general CFT treatment 
in [H]. An analytic derivation of this phenomena remains an open problem. 

5. 4- Spin chain Kondo model 

We now turn to the spin chain Kondo model Eq. (I3.32[) . In previous sections, only the 
uniform part of the impurity part of the von Neumann entropy has been investigated. 
Here we focus on the alternating part which is also present for impurity problems. We 
want to isolate the impurity contribution in the alternating part of S(r,R). Following 
our fundamental definition of Si mp , Eq. (13.34j) . it is also possible to define the alternating 
part of the impurity entanglement entropy: 

Sf mp (J' K ,r,R) = S A (J' K ,r,R) - S A (l,r- 1,R-1). (5.11) 

As before, we have subtracted Sa when the impurity is absent, in which case both r 
and R are reduced by one and the coupling at the end of this reduced chain, linking 
site 2 to 3 and 4, has unit strength. Applying this definition to numerical data involves 
some subtleties. First of all S A is only defined up to an overall sign. Secondly, when 
calculating S^^r — 1,R — 1) we define this as — 5^(1, r', R') with R' — R — 1 since 
the shift from r to r — 1 implies a sign change in the alternating part. For convenience 
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Figure 18. (a) DMRG results for the total entanglement entropy, S(J' K ,r, R) for 
a 102 site spin chain at J$, with a J' K = 0.41 Kondo impurity (■) along with 
S(J' K = 1, r— 1, R—l) (□). For both cases is the extracted alternating part shown along 
with the resulting Sf mp { J' K = 0.41, r, R = 102) for R even, (b) DMRG results for the 
total entanglement entropy, S(J' K , r, R) for a 101 site spin chain at Jf , with a J' K = 0.41 
Kondo impurity (■) along with S(J' K = 1, r — 1,R— 1) (□). For both cases is the 
extracted alternating part shown along with the resulting 5 l ^ lp ( J' K = 0.41, r,R= 101) 
for R odd. For comparison we also show S^, p ( J' K = 0.41, r,R = 102) for i? even from 
panel (a) (dashed line). Reprinted from Ref. [7]. 



we have therefore always exploited this degree of freedom to use a sign convention that 
makes the resulting S^ p positive in all cases. In Fig. [18] we show data for the total 
entanglement entropy along with the extracted alternating parts and the resulting 
for both R = 102 even and R = 101 odd. As was the case for the uniform part of S- imp 
(Fig. El) we do not observe any special features in S^ p (r) for fixed R, J' K associated with 
the length scale and in all cases decays monotonically with r. On the other 
hand, a possible scaling form for was suggested in Ref. [7j. For J' K = 1 we have 
seen above that the alternating part of the entanglement entropy, Sa, is proportional 
to the alternating part in the energy, Ea which, for J2 = J% can be written as follows: 
-EU( r ) = f(j" 7 1 R)l ' \fr f° r some scaling function /. A generalization of the above formula 
to the case J' K 7^ 1 imply that S^^/r should be a scaling function, f(r/R,r/£ K ). 
DMRG results for y/r for fixed r/R are shown in Fig. [19] for a range of J' K and 
R. The values for used to attempt the scaling are the ones previously used for the 
scaling of the uniform part (Fig. [6]). Clearly the results for S^y/r follow the expected 
scaling form. 

6. Impurity entanglement in gapped spin chains 



In previous sections the focus has largely been on work considering critical systems. 
However, impurity entanglement in systems with a gap, corresponding to massive field 
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Figure 19. y/rS^, f° r fixed r/R and a range of couplings J' K at J$. (a) for R 
even and (b) for R odd. The values of £k used are obtained from the scaling of Si mp 
achieved in Fig. [6l DMRG results with m = 256 states. Reprinted from Ref. [7]. 
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theories with a finite correlation length £, have also been considered. In this case the 
generalization of Eq. (11.21) to one dimensional systems with a gap becomes [TT] : 

r 

a / 

So far, mainly two models have been studied, s — 1 chains [BTj 168] at the AKLT 
point [69], E3] an d s — 1/2 J\ — J 2 chains [7J at the Majumdar- Ghosh point, J2 = 
Ji/2 [71]. While the s = 1 spin chain at the AKLT point has a correlation length 
£ = 1/ ln(3), the spin correlations in the s = 1/2 chain at the MG point do not extend 
beyond nearest neighbor and the correlation length is effectively zero. 



6.1. The s = 1 ALKT chain 

Boundary effects in the entanglement entropy of a s — 1 chain was studied in [68] 
building on earlier work [67J. The model considered was the antiferromagnetic s = 1 
Heisenberg chain including a bi-quadratic term: 

L+jV r -l 



H 



E 



-m+i 



Sj ■ S j+ i 



-(Sj ■ S j+1 ) 2 



1. 



(6.2) 



For this special value of the bi-quadratic coupling, termed the AKLT point, the ground- 
state is known exactly [SHI ED] and this fact was exploited to perform exact calculations 
of the entanglement entropy. The entanglement of a sub systems consisting of the 
central section of the chain, 1 . . . L with the remaining left and right parts of the chain 
was considered. Here iVj and N r describe the size of the left and right parts of the 
system and it was shown [68J that the boundary effects in the entanglement entropy 
decay exponentially fast with Ni,N r on a length scale equal to the correlation length 



Entanglement entropy in quantum impurity systems and systems with boundaries 32 



£ = l/ln(3). Hence, impurity effects in the entanglement impurity are in this case 
rather minor. In contrast, by considering a different subsystem impurity entanglement 
in the Majumdar Ghosh model can be rather pronounced. 

The partial concurrence of the two effective s = 1/2 spins at the end of an open 
s — 1 chain has also been studied [30J. 



6.2. The Majumdar Ghosh Model 



B 



(a) 



(b) 



Figure 20. MG state divided in two regions A and B. (a) The total number of site 
R is even and the ground-state is a singlet, (b) \ . Reprinted from [7] 



At the special point J 2 = J\/2, often referred to as the Majumdar-Ghosh [72l[73| [7Tj 
(MG) point, the spin chain model Eq. (13.321) with J' K = 1 is exactly solvable. The spin 
chain has a gap and for R even, in the presence of periodic boundary conditions, a 
two-fold degenerate gound-state of nearest neighbor dimers either between sites 2n + 1 
and 2n + 2 or 2n + 2 and 2n + 3, n > 0, with energy E = —3RJ/8. With open boundary 
conditions, for R even, the ground-state is non-degenerate with dimers between sites 
2n + 1 and 2n + 2, n = . . . (R — l)/2, and with the same energy as the periodic case, 
E = —3RJ/8. See Fig. [201(a). If one now instead considers R odd, an exact form for the 
ground-state wave-function and energy is not known but a very precise variational form 
can be developed [7H E51 [76], E3 EEl El EH] . For R odd it is natural to consider states of 
the following form: 

n 

\n) = \^~^ T - ... -). (6.3) 

Here, — indicates a singlet between site r and r+1 and n — 0, . . . Nd— 1 = (R—l)/2 refers 
to the number of dimers to the left of the soliton, with Nd the total number of dimers. See 
Fig. [201(b). Such states are often called thin soliton states (TS-states) because the soliton 
resides on a single site and is not 'spread' out over several sites as would have been the 
case if one had included states with valence bonds longer than between nearest neighbor 
sites. Note that the soliton only resides on the odd sites of the lattice, r = 1, 3, . . . , R. It 
is important to realize that the TS-state as defined are not orthornormal and even though 
it is straight forward to form linear combinations \m) that are orthogonal [7HE9] through 
the transformation \m = 0) = \n — 0), \m) = (2/y/3) [\n) + (1/2) \n — 1)] m > 1, one 
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finds that for the purpose of calculating the entanglement this orthogonalization is less 
useful at the initial stage of calculating the entanglement entropy. 

If only the ^ state of the gound-state doublet for R odd is considered one can write 
a thin soliton ansatz (TS-ansatz) for the ground-state wavefunction: 

N d 

l4s>^ECV)- (6-4) 

n=0 

It can be shown that the restriction to the " state is not important since any 
linear combination of the degenerate ^ and ^ states will yield the same entanglement 
entropy [TJ. The components of the ground-state wavefunction, ip^ 1 can be obtained in 
a variational manner or through a simple analytical estimate [7J : 



C°Wi^sin(^±S). (6.5) 



N d + 2 \ R + 3 

With the i\) s ° l determined we can proceed with a calculation of the entanglement entropy. 
Due to the dimerization of the ground-state one sees that for R even (and open boundary 
conditions) S(r, R) simply oscilates with r between ln(2) and depending on whether r 
is a site at the beginning or end of a dimer, respectively. With a J' K ^ 1 in Eq. (I3.32p the 
entanglement is much richer and exact results are not available, however, very precise 
variational calculations [7J of S(r, R) are possible. Let us first consider R odd and 
J' K = 1 with r odd. We first divide the system in two parts (A, B) at r and let \<p) 
denote a basis for 1 . . . r and a basis for r + 1..R. We then have: 

3 

I*ts) = ^CuWiMj). (6-6) 

i,j=0 

The above form follows since the system was divided in two parts at r and since r is 
odd. If the soliton is to the left of r and since r is odd the division between A and B 
will not 'cut' a dimer. However, if the soliton is to the right of r the division will cut a 
dimer and we effectively get an additional soliton at the end of the A space leading to 
the following 3 separate cases for 



hM = I>n1^~^ T ■•••-) 
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iik> = r- - - ... — i). (e.?) 

With the corresponding definitions for \4>j): 
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(a) (b) 

Figure 21. (a)DMRG results withm = 256 for S(J' K = l,r,R = 201) for the 
spin-chain model with J2 = 1/2 (MG model). The small solid circles represent the 
theoretical result obtained using the TS-ansatz. Also shown arc the uniform parts 
Su{J' K = l,r,R = 201) and S V {J' K = l,r,R = 200) the difference of which is 
S imp (J' K = 1). (b)DMRG results withm = 256 for S(J' K = 0,r,R = 200) for the 
spin-chain model with J 2 = 1/2 (MG model). The small solid circles represent the 
theoretical result obtained using the TS-ansatz. Also shown arc the uniform parts 
Su(J' K = 0,r,R = 200) and S V {J' K = l,r,R = 199) the difference of which is 

<Simp(^R" = 0). 



— 2 

I0 3 ) = £ rA +n \ T T - ... -)• (6.8) 

n=0 

With these definitions it immediately follows that: 

l*$S> = llMIM + 4g Dlfe>l&> - 1^3)103)] • (6.9) 

from which the coefficients C defined in Eq. (16. 6ft can be determined. The states \(f>i) and 

are clearly not orthonormal, however, by explicitly orthonormalizing the states 
region B can be traced out and the reduced density matrix p for region A determined It 
is straight forward to generalize this approach to the case where r is even (and R is odd 
with J' K = 1). Results from such a calculation is shown in Fig. l2lTa) where they are 
compared to DMRG results. Excellent agreement is observed. For comparison, results 
for S U {J' K = l,r,R = 201) and Su{J' K = l,r,R = 200) are also shown in Fig. EH(a). 
The difference of these two uniform parts yields Simp- The influence of the impurity 
spin coupled with J' K = 1 is clearly visible and extends over the entire range of r. 

As already outlined in section 13.21 the entanglement of the impurity spin with the 
bulk of the chain is sizable (one might even say maximal) even when J' K = 0. With 
some additional algebra it is possible to extend the thin-soliton approach also to this 
case [7] by considering a decoupled impurity entangled with a bulk chain of odd length 
R — 1. Hence, the TS-ansatz can be applied to the bulk chain and precise variational 
results obtained. It is crucial to note that the impurity spin and the bulk chain form 
a singlet state even though their coupling, J' K , is zero. That is, we're considering the 
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J' K — > limit of the entanglement entropy. Results obtained from the TS-ansatz are 
shown in Fig. [21Tb) where they are compared to DMRG results obtained using spin- 
inversion. The DMRG calculations are performed on systems where the impurity spin 
is explicitly included resulting in a 4-fold degenerate ground-state consisting of a singlet 
and a triplet. The application of spin- inversion symmetry is therefore necessary in 
order to distinguish the singlet state of interest from the degenerate triplet state. The 
agreement between the numerical and variational results in Fig. l2~iT b) is excellent, the 
small discrepancies visible at a few values of r are due to complications using spin- 
inversion [80] in the DMRG calculations specific to this value of J 2- For comparison, 
results for S U (J' K = 0,r,R = 200) and Su(J'k = l,r,R = 199) are also shown in 
Fig. l2TT b). The difference of these two uniform parts yields S[ mp which is non-zero 
throughout the system. 

Perhaps surprisingly, it is possible to obtain a much more intuitive picture and yet 
still very precise by simply assuming that the states and are orthogonal. With 
this assumption it is easy to see from Eq. (I6.9P that the reduced density matrix for 
region A (with R odd) is simply: 
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Here p = X^n=o l^/C'l 2 ~ rjR — sin(27rr/i?)/(27r) is the probability of finding the 
soliton in region A. From these expressions for p the entanglement entropy can easily 
be evaluated. In order to calculate 5i m p the uniform part of 5" for an even length chain 
with J' K = 1 is needed, but, as mentioned previously, this is simply ln(2)/2 and it follows 
that: 

S imp (J' K = hr,R) = -pHp) -(1-p) ln(l -p) = Sf™, R odd. (6.11) 

In this case the impurity entanglement arises solely from the entanglement of a single 
particle (the soliton) that is present in the ground-state and one therefore refers to this 
contribution as the single particle entanglement, Sf^F. See Fig. [22T b). 

With more effort an analoguous calculation can be carried through for the case of 
R even and J' K = 0. One finds [7]: 

SimM = 0, r, R) = {1- p) ln(2) = ^ p B , R even. (6.12) 

We see that in this case there is no contribution from the single particle entanglement 
as one would expect since there are no solitons (single particles) present in the ground- 
state. Instead, the impurity entanglement is given purely by the impurity valence bond 
impurity valence bond (IVB) (See section 13. 21) where in the present case one identifies 
the probability that the IVB does not cross the boundary between regions A and B 
with the probability of finding a soliton in region A. See Fig. [22T a). 

Some surprising observations can be found by performing numerical calculations 
of "S'imp away from the MG point where the above variational calculations are quasi- 
exact. One finds [7J: (i) S^mp for both even and odd R is non-zero over the entire 
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Figure 22. (a)DMRG results withm = 256 for S(J' K = 0,r,R = 400) for the 
spin-chain model with J2/J = 0.2, J|, 0.3, 0.4, 0.5 shown along with the result for 
Simp = (1 -p)ln( 2 ) with P calculated at J 2 /J = 0.5 (dashed line). (b)DMRG 
results withm = 256 for S(J' K — l,r,R — 401) for the spin-chain model with 
J 2 /J = J|, 0.3, 0.4, 0.5 . The dashed line represents Sf™ = -pln(p) - (1 -p) ln(l-p) 
with p calculated at J2 = J/2. Reprinted from [7]. 



range of r and is not limited by the correlation length which is effectively zero, (ii) 
Sim P (J' K = 0,r,R) for R even changes only slightly when J 2 is decreased from J 2 = J/2 
to J 2 = Jl and for J 2 < J| in the gap less phase it appears not to change at all with J 2 . 
In all cases, < J 2 < J/2, does the IVB picture seem to correctly describe S^. This is 
illustrated in Fig. 1227 a). IT is perhaps surprising that the IVB picture works relatively 
well for J 2 < J| where long-range valence bonds are present in the ground-state, (iii) 
S\xa V {J' K = l,r,R) ~ Si™ f° r ^ °dd decreases rapidly with J 2 from J/2 to J2 where 
it vanishes. For J 2 < J| this part of the impurity entanglement is negligible. Hence, 
it seems likely that this fact is related to the system becoming gapless at Jf. This is 
illustrated in Fig. l22Tb). 

Finally we note that, the concurrence of the end spins in the dimerized J\ — J 2 
s — 1/2 model have been studied [30J. 

7. Conclusions 

We have reviewed recent results related to the impurity contribution to the 
entanglement, Si mp , arising from a quantum impurity or boundary. Most notably it 
has by now clearly been established that the entanglement can be dramatically changed 
by the presence of an impurity even in the case where the physical coupling to the 
impurity is zero. The role played by different boundary conditions in 1 + 1 dimensional 
critical systems is well understood and in agreement with results from CFT. However, 
for the case of mixed boundary conditions (section 13. 2j) quantitative theory is not yet 
available and a detailed understanding of the R dependence of bulk impurity effects 
(section HJ) would be desirable. A consistent picture of the entanglement of a qubit 
with the environment as described by the spin-boson problem and the entanglement 
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arising from Kondo impurities has been developed based on established theory of Kondo 
systems. A more heuristic picture of impurity entanglement in gapped (and to a certain 
extent also critical one dimensional systems) based on valence bond physics and matrix 
product states is emerging. Useful notions of single particle entanglement (SPE) and 
impurity valence bonds (IVB) have been introduced. Some details of this picture are still 
missing. As an example, how the single particle entanglement disappears as the system 
approaches criticality is still an open problem. Comparatively few results are available 
for quantum impourity entanglement in two (and higher) dimensional quantum systems 
and a detailed theory is still lacking. 
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